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Abstract 

j^ '. The decay widths of S-wave quarkonia {r]c, r]b — > 77 and J/V', T" 



in 



(N 



e"'"e~) are calculated on the basis of a quasipotential approach. The 
nontrivial dependence on relative quark motion of decay amplitude is 
taken into consideration via quarkonium wave function. It is shown 



^ ' that relativistic corrections may be large (10-50 %) and comparable 

O 



CS ■ with QCD corrections. 

o 



1 Introduction 

CT^ I The constituent quark models O are being not well grounded from the point of view of 

quantum field theory, as the QCD sum rules method 0] or lattice QCD approach 0], but 

C^hI they give us a opportunity to describe most of all existing experimental data. The great 

^ success in description of heavy meson static characteristics and decay widths has been 

g_) I obtained on the basis of a relativistic quark models Q, ^. As it was shown, relativistic 

,S^ • corrections are not small for mesons containing comparably light charm quark f^. In 

this paper we present the results of calculation of the relativistic corrections in S-wave 

. ■ quarkonium decays {ric,r]b — * 77 and J/ip,T —^ e'^e^) using the approach based on a 

rS I Logunov-Tavkhelidse local quasipotential equation [Q and pcrturbative QCD methods. 



2 Decay amplitude in quasipotential approach 

The meson decay amplitude m.ay be presented as the product of an parton decay ampUtude 
Mh{qiiQ2) with amputate quark hues times a Bethe-Salpeter wave function ^(gi, (72): 



M 






Sp[^{qi,q2)MH{qi,q2) 



(1) 



where ^{qi, (72) describes two-quark bound state, A^H(gi, 92) describes transition between 
initial two-quark state and final state of free particles, Sp means the trace over spin and 
colour indexes. The wave function '^{qi, 92) satisfies two-particle Bethe-Salpeter equation, 
which has not physically interesting exact solutions. That is why the redefinition of the 
wave function for calculation of decay amplitude A1 is needed. We must to introduce new 
wave function, which has clear physical interpretation and it is determined by more simple 
bound state equation. 

Let us define vertex function T(qi, (72) as follows: 



*(gi,92) = — ^—T{qi,q2)— , , . . 



(2) 



Transforming (1) to the center-of-mass reference frame and using projection operator 
decomposition for quark propagator on positive and negative energy states 



p + m 



1 



U°'{p)U°'{p) V'^{-p)V°'{-p) 



p^ — m? + 10 2e(p) \_p° — e{p) + ; 
we can write amplitude A4 in following form: 
1 f (Pq f dq° 



M 



p° + e{p) — io 
{M-2e{q)) 



(3) 



(27r)3 7 2e{q) J 2Tr {q° - e{q) + io){q° - M + e{q) - io) 
Sp\'i'Miq)Rsiq)MHiqi,q2) 



(4) 



where e{p) = ^^fP + w? , qi = {q°, q), (72 = {M — q°, —q), M is the mass of bound state, 
^A/(<?) is single time bound state wave function which was obtained via projection the 
Bethe-Salpeter amplitude on positive energy state [g[: 

Uiq)T{qi,q2)V{q) = 2e{q){M ~ 2e{q))^Miq). (5) 

Note, that for S-wave bound states which are discussed here, 'i'M{q) — ^J\/(|'7|)- 



The spin property of the two-quark bound state is described by relativistic projection 
operator, in which the exact dependence on vector momentum g must be considered. We 
write it as in H: 

j^ (g + m)(l+7o)g(m^gO 

^^' 2V2{e{q)+rn) ^^ 

where s = e — 6^7'^ for vector ^Si state, e^ is polarization vector, and s = 75 for 
pseudoscalar ^So state, q = {e{q),q), q' = {e{q),—q). 

The relative energy q^ integral in (4) is defined by poles of integrated function. The 
pole q° ^ e{(j) — iO in lower halfplane of complex variable q^ gives main contribution. It 
is obviously that nontrivial contribution of other poles in parton amplitude A^_f/(qi, 92) is 
suppressed by factor A « {M — 2e{q))/m, which however is not very small in the case of 
charmonium state decays and it will be discussed separately. 

The nonrelativistic approximation corresponds to limit (f — > 0. In this case the meson 
decay amplitude is equal to 



Mo = ^— Sp 



RsoMh{0,0) , (7) 



where 

*m(0) = f T^^^^Af (g), and Rm = m(l + 7o)s/%/2. 

The expression (7) is agree with usual used formula in potential quark models m. 

The wave function ^m(p) of two-particle bound system with mass M satisfies to 
quasipotential equation in the center-of-mass reference frame: 

(M - 2e{p))^M {p) = J 703 V{p, q, M)^M (q) , (8) 

where V{p, q, M) is so-called quasipotential, which is calculated using two-particle scatter- 
ing off-shell amplitude, p is a relative momentum of particles: p = (pi — p2)/2. For more 
comfortable investigation bound state wave function and mass spectrum, we transform eq. 
(8) to a local form [|lO): 

b\M) f \ ^^^^^ ^ I _A_^(^. . M)^M{q), (9) 



2y.R 2^uJ J (27r) 



where b^{M) — {Iv'P — Am?)/ A is the square of relative momentum of particles on energy 
shell (M — 2e{p)), /i/j = M/4 is relativistic reduced mass. Than we obtain from eq.(9) 
the bound state equation in coordinate space: 

^ +[/(f))*Af(r=) = ^fl*Af(r), (10) 



where 



bl NP - An? 

R 



E--^ 



2fiR 2M 

In the case of spherical symmetry potential U = U{r) we have for radial part of wave 
function: 

f2^-2f^n{U{r)-En)x^0, (11) 

where 

x{r) — vATTr'^Mir). 

Taking into account the accuracy of our calculations, we shall use non-relativistic 
approximation of quark-antiquark interaction potential. As usual we present it using a 
composition of Coulomb and linear components: 

[/(r) = har + c. (12) 

r 

In our calculation we use following numerical value for parameters of interactive po- 
tential: b = 0.25, a = 0.27 GeV^ c = -0.76 GeV - for c-quarks; 6 = 0.4, a = 
0.25 GeV^, c = —0.3 GeV - for 6— quarks ^. Because we don't take into considera- 
tion spin dependence of quark-antiquark potential, it is necessary to use quark masses as 
free parameters for right description of mass spectra of quarkonium states. We obtained: 
TOc = 1-62 GeV and rub = 4.87 for ^Si state, m^ = 1.55 GeV m^ = 4.84 GeV for 
1^0 state. In this case we have: Mj/^, ^ 3.10 GeV, M^^ = 2.98 GeV, Mr = 9.45 GeV, 
Mr,, = 9.4 GeV. 

So, solving numerically bound state equation (11) with potential (12), we obtain 
quarkonium wave function in coordinate space. After numerical Fourier transformation 
we obtained the wave function in momentum space, which was used for estimation of 
relativistic effects in S-wave quarkonium decays. 



3 Two-lepton decays J /ip^ T — > e^e 

Two-lepton decay of '^ Si state of heavy quark- antiquark bound system ( for example J /iji) 
is very interesting both for theory and experiment in particle physics. First, two-lepton 
mode of J/V' decay is good trigger in experiment, where J/V' production is studied. Sec- 
ond, as it is usually expected, we can used in first approximation simple non-relativistic 
description for bound state of heavy quarks and study relativistic effects as small cor- 
rections. From the point of view perturbative QCD and motivation of quark-antiquark 
potential, the small value of running constant q;s(to^) = 0.3 is also very important. That 
is why calculation J/tj) or T decay widths is a good test of our knowledge about quark- 
antiquark interaction at large distance as well as about role of relativism in description of 
heavy quark bound states. 

The amplitude of decays J/^/", T — > e+e~ has the form: 

M{J/^^e+e-) = V2M I ^-£^ J ^Sp[Rjiq)MHiqi,q2)] 

p (M-2£(g-))Vl/M(g) .^3. 

" {q° - e{q) + io){q° - M + e{q) - ioY ^ ' 

where qi = {q°,q), q2 = {q° — M,—q), F^ ~ \/3 is colour factor. The relativistic 
normalization for quasipotcntial wave function gives the multiplier y2M in (16). The 
"hard" part of amplitude can be written as: 



MH{quq2) = ^rUeih)-f^Veik2). 



In the case of J/ijJ, T -^ e+e^ decays, amplitude A4h doesn't depend on relative quark 
energy and the q° integration is carried out simply. After sum over lepton polarizations, 
the squared modulus of amplitude takes the form 



(27r)'^ e{p){e{p) + m 



2 _ e^el f <fq *m(<L) 



M3 J {2T:f e{q){e{q) + m) 

f (tX r->f7tl , c,Aki,k2)n^iq)n-'ip), (14) 

J (2TTy s{p}(e{p) + m) 



where 



Cf^„{ki, fc2) = 4 [fci^fc2^ + ki^k2fj^ - (fcifc2).9pi.] , 
H'^iq) = Sp [(<7i + to)(1 + 7o)e,/('Z2 - m)^] ■ 



Using the condition for J /tp polarization vector 



we average (14) over J/^p polarizations and than carry out the angle integration. Finally 
we obtain: 

16Q=^eg 
3^^'" "' 9n^M 



where 



/n = 



qx{q)dq 



(m^/o + 3m^/i + 2/2) 



h = 



qx{q)dq 



,h^ 



£{q)qx{q)dq 



£{q){£{q) + -m)' J {e{q)+m)' J {e{q) + m) 

The width of Jf-ip -^ e+e^ decay and amplitude are connected as follows: 



T{J/ijj -^ e+e-) = 



' \T.\M\' 



WttM 3 



In non-relativistic approach it reads: 



2 2I ^m(0)P 



(15) 
(16) 

(17) 
(18) 



Our numerical calculation gives for J /^ particle F/Fo — 0.89 and for T particle - 
F/Fo = 0.97. So, the values of relativistic corrections in two-lepton decays of J/V', T" 
particles are not small, but they are less than QCD corrections in next order of a^ [O]: 



^QCD — To 1 — 



16 Os 



(19) 



Putting up as{m1) — 0.3 and as(TO^) = 0.2, we obtain TQC'D/^o{J/fp) ~ 0.5 and 
rQC£>/Fo(T) w 0.7. Note, that our calculations give more large relativistic corrections 
than estimation in binding energy e approximation H : 



F 



1 e 
3 m 



(20) 



This equation gives F/Fq = 0.96 for J/'ip and F/Fq — 0.98 for T particles. The difference 
may be explained using following fact. We exactly take into consideration relative quark 
motion in amplitude and in wave function, but in B fp' /m? approximation in amplitude 
is used and wave function is presented only in origin. 
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4 Two-photon decays ?7c, r]b 



77 



The two-photon decays of rjc and rji, particles are not studied sufficiently in experiment 
opposite to lepton decays of J/ip and T particles. The experimental uncertainty of decay 
width is equal to 50% for rj^ ~^ "f"f decay and data is absent for 77;, -^ 77 decay p2| . 
Contradict to this fact, these decays are very interesting for discrimination of many models 
which describe heavy quark bound state and QCD properties at large distance [ |l3| . 
The decay amplitude rjc -^ 77 {r/i, — > 77) may be presented in the form: 



M{r]c^jj) = 



^2M 



d\ 



^ J ^Sp[Rr,{q)MH{qi,q2)] 



Fr 



2e{q){2 

{M-2e{q))^M{q) 



■ {q" - e{q) + io){q° ~ M + e{q) - io) ' 
The "hard" part of amplitude Ai{ijc -^ 77) may be written as follows: 



(21) 



MH{qi,q2) = e^e^e^(A:i)e^(A:2)[7'' 
qi-k2-\-Tn 



qi - fci 



{qi - h) 



,r + 



{qi - k2Y 



(22) 



The integral function in (32) has two poles in plane of complex variable q°: q^ — e{q)—io 
or q° = M — e{q) + io as well as (7° = fc° ± y^k^"^ + m? ^io and g^' = ^2 ^ \/kf^ + rn^ T io . 

The contribution of other poles from " hard" part of amplitude is suppressed compa- 
rably the contribution of pole g" = e{q) by factor A w 0.1 for r/c and A w 0.03 for rji, 
particles. We will discuss it below. 

After integration equation (32) over gj', the contribution of the pole q^ — e{q) — io 
reads: 



M{iic -^ 77) = V2M 



Sq ^MJq) 
(2^)3 2e{q) 



F,Sp R^{q)MH{q) 



\M\'^i2F^e^,{A7Ta) 



,MJ^ 



where 



J 



(27r)3 



d^q T .^ 1 , 

2e{q) q 



Am? 
e{q) 



(23) 
(24) 



e{q) 



In non-relativistic approach we obtain well known result from (32): 

|A4„P = 32F,M(4.af^:^. (25) 

After numerical calculation of (35) we obtained T/Tq = 0.53 for rjc -^ 77 decay and 
r/Fo = 0.71 for rih — > 77 decay. The increase of the relativistic corrections in two-photon 
decays as compare with two-lepton decay is explained by the nontrivial dependence of the 
two-photon decay amplitude from relative quark momentum q. Obtained value of rela- 
tivistic correction for two-photon decay width is approximately equal to QCD correction 

/ as, 20 7r2 
rQci5-r<,(^l --(---) 

For estimation of contributions from the other poles in "hard" part of decay ampli- 
tude, we keep in mind that we can to ignore the dependence of amplitude on the vector 
momentum q. Thus the given contribution elementary related with non-relativistic ap- 
proximation: 

M2=Mo{l + S), (26) 

where 

, M - 2m 

6 = — -= . 

v8to 

For rjc — *■ 77 decay we obtained: S — —0.027. Note, that the sign of this correction is 
strongly dependent on choice of quark-antiquark potential. 

5 Conclusions 

Using relativistic quasipotential approach, we have calculated decay widths for S-wave 
quarkonia: rjcrji, — > 77, J/^,T — > e~^e~. Nontrivial dependence on relative quark 
motion of decay amplitude is taken into account both in parton amplitude as in quarkonium 
wave function. 

Our conclusions arc: 

1. Relativistic corrections for S-wave quarkonium decay are not small (10-50%) and 
comparable with QCD corrections. 



2. Quasipotential approach gives more large value for relativistic corrections than bind- 
ing energy approximation. 

3. Additional poles in a relative-energy integral determine non vanishing contribution 
in decay width of quarkonia. 
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